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Abstract 

In this article a mathematical model of an immobilized 
oxidase enzyme electrode is presented. The model is based 
on the three reaction-diffusion equations containing a non- 
linear reaction term under the steady state conditions. A 
simple analytical expressions pertaining to concentrations of 
the immobilization of three enzyme substrates are obtained 
by using the new Homotopy perturbation method (HPM). A 
simple analytical expression of the concentrations of 
substrate, oxygen and oxidized mediator and current was 
obtained in terms of the thiele moduli and the small values 
of the normalized surface concentration of substrate B s , 
thenormalized surface concentration of oxygen B 0 and the 
normalized surface concentration of oxidized mediator B m . 
These analytical solutions are compared with the numerical 
simulation. A good agreement between analytical 
expressions and numerical results is noted. 
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Introduction 

The basic concepts of enzyme electrodeswereintrodu- 
ced by Clark and Lyons [Clark et. al (1962), Lyons et.al 
(1962)]. Many kinds of biosensors based on 
electrochemical enzyme electrodes have been studied 
and a majority of this device operates in an 
amperometric mode [Scheller et.al (1992), schulbert et. 
al (1992)]. Biosensors are usually classified into various 
groups either by type of transducer employed 
(electrochemical, optical, piezoelectric, and thermal) or 
by the kind of bio-recognition element utilized 
(antibody, enzymes, nucleic acids, and whole cells). 
Both components of the biosensor, namely, the bio 


recognition element (referred as a receptor) and 
transduction platform (referred as a transducer) play 
an important role in the construction of a sensitive and 
specific device for the analyte of interest. 

Schulmeister has described models for multilayerand 
multienzyme electrodes; these models assumed 
operationof the electrode under diffusion control, such 
that the enzymekinetics are linear with substrate 
[Schulmeister et. al (1990) and Pfeiffer et. al (1993)]. 
The kinetics is described by a parabolic differential 
equation with linear inhomogeneities. Thefact that 
relatively few enzyme electrodes have been 
commercialized may be due to the technical 
problemsassociated with either the enzyme reaction or 
the underlying product/substrate-sensitiveelectrode, 
such as the availability of appropriate enzymes, their 
successful immobilization, andprobably most 
importantly the achievement of an acceptablecatalytic 
lifetime of the enzyme. 

Leypoldt and Gough have described a model for a two 
substrate enzyme electrode with the non-linear 
enzyme reactionfLeypoldt et. al (1984) , Gough et.al 
(1984)]. This model to describe the behaviorof a glucose 
oxidase electrode. Bergel and Comtat describe the 
transient response of a mediated amperometric 
enzyme electrode by using an implicit finite difference 
method [Bergal et. al (1984), Comtat et. al (1984)]. A 
two-substrate enzyme electrode was described by 
Gooding and Hall [Gooding et. al (1996) , Hall et. al 
(2012)]. The development of models for enzyme 
electrodesprovides a better understanding of the 
individual processesinfluencing the response of the 
device, and this informationmay be used as a guide for 
directions for improvement of thesensor design. 
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The importance of enzyme electrode and the various 
approaches to enzyme electrochemistry has been 
discussed. The current response has been also 
observed experimentally for mediated enzyme 
electrodes employing glucose oxidase [Martens et. al 
(1995), Hindle et. al (1995), Hall et. al (1995) and 
Pallachi et. al (1990), Turner et.al (1990)]. The 
membranes provide an ideal support for the 
immobilization of the biocatalyst. Substrate partition 
at the membrane/fluid inter-phase can be used to 
improve the selectivity of the catalytic reaction 
towards the desired products [Trevan et. al (1981)]. 
Homogeneous membranes are used as carries for 
immobilization of enzymes [Simon et. al (2002), 
Hailweil et. al (2002), Cass et. al (2002)] and also used 
in biomaterials, bio-separators and biosensors [Liu 
et.al (1996) ,Zhang et. al (1996) , Yu et.al , Deng et. al 
(1996) ]. 

In this paper, we describe a model for thisthree- 
substrate enzyme electrode and employ the model 
toinvestigate the influence of the oxygen on the 
current responseof a mediated oxidase electrode. To 
our knowledge, no general analytical expressions that 
describe concentration of substrate, oxygen and 
oxidized mediator for various of the thiele modulus cp 2 , 
cp 2 and <p 2 m and the normalized parameters B s , B 0 and 
B m have been reported. However, in general, 
analytical solutions of non-linear differential equations 
are more important manipulation and analysis. For 
this reason, we have derived that analytical 
expressions corresponding to the concentrations of 
substrate, oxygen and oxidized mediator in an oxidase 
enzyme electrode using new Homotopy perturbation 
method. This method is reliable and highly accurate in 
handling non-linear problems. 

Mathematical Formulation of the Boundary 
Value Problem 

The details of the model adopted have been fully 
described in Martens and Hall [Martens et. al (1994), 
Hall et.al (1994)]. Figured represents the general 
kinetic reaction scheme of enzyme-membrane 
electrode geometry [Gooding et. al (1996), Hall et.al 
(1996)]. The general reaction scheme for an 

immobilized oxidase inthe presence of two oxidants 
can be written as follows: 

Eax+S-L+ES-^E^+P 


E red + ^2 ' * ^ OX + ^2 ^2 


(1) 

(2) 


E red + Med ox — E ox + Med red (3) 


We assume that the concentrations of all reactants and 
enzyme intermediates remain constant for all time. 
Also the concentration of total active enzyme \E t ] and 
the reactants in the bulk electrode remain constant. We 
can consider that the diffusion of the reactants can be 
described by Fick's second law and the enzymes are 
assumed to be uniformly dispersed throughout the 
matrix. The enzyme activity is not a function of 
position. The reaction/dif fusion equations correspond- 
ing to the concentrations of substrate, oxygen and 
oxidized mediator within a matrixcan be expressed 
as[14] 
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where [S ] , \Med 0X \ and [0 2 ] are the concentrations 

of substrate, oxidized mediator and oxygen and y is 
distance from the electrode. D s , D m and D 0 are the 
diffusion coefficients of substrate, oxidized mediator 
and oxygen. [E t ] Denotes the concentration of the 
total active enzyme in the matrix and J3 S - [k 2 + k 3 ) / k x , 
P Q = k 3 / k 4 & fl m =k 3 / k 5 represent reaction constants 
of substrate, oxygen and oxidized mediator 
respectively. The boundary conditions are given by 

When y -d , 

[ 02 \ = [0 2 ] b =K 0 [0 2 l,[S] = [S] b =K s [S] b , 
[Med ox ] = [Med ox ] h = K m [ Med ox ] x (7) 

when y = 0, 

[Med 0X ] = [Med 0X \ d[0 2 \/ dy = d[S]/ dy = 0 (8) 

Here [S\, [0 2 \ [Med ox ] b and[S„] , [0 2 l , [Med ox ]^ 

are the bulk concentrations of substrate, oxygen and 
oxidized mediator species. K 0 , K s and K m are the 
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partition coefficients for oxygen, substrate and 
oxidized mediator respectively. Equation (4) can be 
represented in the normalized form using the 
following normalized parameters. 


_[ s ] ^[<>2] r _[ Med ox] _y 

"[*]/ °~[QiV m ~[Med ox } b ' X ~d ’ 
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where F s , F 0 and F m represent the normalized 
concentrations of substrate, oxygen and oxidized 
mediator and B s , B 0 and B m are the corresponding 
normalized surface concentrations. The surface 
concentration is the ratio of the bulk concentration and 
the reaction constants. <p 2 , <p 2 and cp] n denote the 
thiele moduli of substrate, oxygen and oxidized 
mediator, respectively. Thiele modulus cp 1 represents 
the ratio of the characteristic time of the enzymatic 
reaction to that of substrate diffusion and dis the 
thickness of the enzyme layer. The systems of three 
non-linear reaction/diffusion equations in normalized 
form are 
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The boundary conditions become 

F=F=F„=latx = \ 


s o m 


dF dF 

F m =1, — - = — — = 0 < 2 / x = 0 
dx dx 

The normalized current J ox is given by. 


J ox ~ 


d_K 
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(15) 

(16) 

(17) 


Analytical Expression of the Normalized 
Surface Concentrations Using the New 
Homotopy Perturbation Method 

Linear and non-linear phenomena are of fundamental 
importance in various fields of science and 
engineering. Most models of real - life problems are 
still very difficult to solve. Therefore, approximate 
analytical solutions such as Fiomotopy perturbation 
method (HPM) [Gohri et. al (2007), Ozis et. al (2007), 
Li et. al (2006), Mousa et. al (2008), He (1999 and 2003), 
Ariel (2010), Ananthaswamy et. al (2012 and 2013)] 
were introduced. This method is the most effective 
and convenient ones for both linear and non-linear 
equations. Perturbation method is based on assuming 
a small parameter. The majority of non-linear 
problems, especially those having strong non-linearity, 
have no small parameters at all and the approximate 
solutions obtained by the perturbation methods, in 
most cases, are valid only for small values of the small 
parameter. Generally, the perturbation solutions are 
uniformly valid as long as a scientific system 
parameter is small. However, we cannot rely fully on 
the approximations, because there is no criterion on 
which the small parameter should exists. Thus, it is 
essential to check the validity of the approximations 
numerically and/or experimentally. To overcome these 
difficulties, HPM have been proposed recently. 

Recently, many authors have applied the Homotopy 
perturbation method (HPM) to solve the non-linear 
boundary value problem in physics and engineering 
sciences [Ghori et. al (2007), Ozis et. al (2007), Mousa 
et. al (2008)]. Recently this method is also used to solve 
some of the non-linear problem in physical sciences 
[He 1999 and 2003 ]. This method is a combination of 
Homotopy in topology and classic perturbation 
techniques. Ji-Huan He used to solve the Lighthill 
equation [Mousa et. al (2008)], the Diffusion equation 
[He 1999] and the Blasius equation [He 2003]. The 
HPM is unique in its applicability, accuracy and 
efficiency. The HPM uses the imbedding parameter p 
as a small parameter, and only a few iterations are 
needed to search for an asymptotic solution. The 
simple analytical expressions of the concentrations of 
substrate, oxygen and oxidized mediator are as 
follows: 


^W = y6 2 -i) + 1 

(18) 

^W=fb 2 - l ) +1 

(19) 


34 



Frontiers in Sensors (FS) Volume 2, 2014 


www . seipub . or g/f s 


where. 
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y {B 0 + B m )B s+ B s+ (B 0 + B m )j 

The corresponding normalized current using eqn. (20): 


/ - — 
Jox 2 


(24) 


where k\ is defined by the eqn. (23). 


increases, the corresponding normalized concentration 
of substrate F s , oxygen F 0 and mediator F m decreases 

for some fixed values of B 0 , B m when $ 2 S =1, <p* = 1 
and <p 2 m = 1 . 

Figure4 is the normalized concentrations of 
(a)substrate F s , (b)oxygen F 0 and (c) mediator F m 
versus the normalized distance x. From Fig. 4 (a)-(b) it 
is clear that when the normalized parameter B 0 
increases, the corresponding normalized concentrations 
of substrate F s , oxygen F 0 decreases for some fixed 

values of B s , B m when^ 2 s =l, (p 2 0 - 1 . Figure. 4 (c) it is 
clear that when the normalized parameter B s 
increases, the corresponding normalized concentrations 
of mediator F m increases for some fixed values of B s , 

B m and q>l= 1 . 


Numerical Simulation 

The non-linear reaction-diffusion equations eqns. (12)- 
(14) for the boundary conditions eqns. (15) and (16) are 
also solved numerically. We have used the function 
pdex4 in Scilab/Matlab numerical software to solve 
numerically, the initial-boundary value problems for 
parabolic-elliptic partial differential equations. This 
numerical solution is compared with our simple 
analytical equations in Figs. (2) - (6). A satisfactory 
agreement is noted for various values of the thiele 
moduli cp 2 s , (pi and (fl and possible small values of 
the dimensionless parameters B s , B 0 and B m . 

Results and Discussions 

Figured is the schematic model of an enzyme- 
membrane electrode. Figure2 is the normalized 
concentration of (a)substrate F s , normalized 
concentration of (b) oxygen F 0 and normalized 
concentration of (c) mediator F m versus the 
normalized distance x. From Fig.2 (a), (b) and (c) it is 
clear that when the thiele modulus cp 2 s , cp 2 0 and (p 2 m 
increases, the corresponding normalized 
concentrations of substrate F s , oxygen F 0 , mediator 
F m decreases for some fixed values of the normalized 
surface concentrations of the substrate B s , normalized 
surface concentration of oxygen B 0 and normalized 
surface concentration of the mediator B m respectively. 

Figure3 is the normalized concentrations of 
(a)substrate F s , (b)oxygen F 0 and (c) mediator F m 
versus the normalized distance x. From Fig. 3 ((a)-(c)) 
we infer that when the normalized parameter B s 


Figure. 5 is the normalized concentrations of 
(a) substrate F s , (b) oxygen F 0 and (c) mediator F m 
versus the normalized distance x. From Figs. 5 (a) and 
(c) it is clear that when the normalized parameter B m 
increases, the corresponding normalized concentrations 
of substrate F s and mediator F m decreases for some 

fixed values ofi^, £ 0 when (j) 2 s - 1 and (p 2 m =1. Figure. 
5 (b) it is clear that when the normalized parameter B m 
increases, the corresponding normalized 
concentrations of oxygen F 0 increases for some fixed 

values of B s , B 0 and cp 2 0 = 1 . 

Figure. 6 is the normalized concentrations of the 
substrate F s , oxygen F 0 and the mediator F m versus the 
dimensionless distance x . From this figure, we note 
that the substrate, oxygen and mediator increasesfor 
some fixed values of the parameters B s ,B 0 ,B m and 

cp 2 s = <pl = (pi = 1 Figure. 7 is the variation of normalized 
current J ox with (a) normalized thiele modulus (p 2 m for 
the mediator (b) normalized surface concentration of 
the substrate B s (c) normalized surface concentration 
of oxygen B 0 and (d) normalized surface concentration 
of mediator B m . From these Figs. 7 (a), (c) and (d) the 
normalized current decreases and (b) increases for 
some fixed values. 

Far wall Membrane Electrode Solution 



d » 


FIGURE 1: GEOMETRY OF THE PROBLEM 
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(a) 



(b) 



(c) 



FIGURE 2: NORMALIZED CONCENTRATIONS OF (A) 
SUBSTRATE^ (EQN. (18)) (B) OXYGEN F 0 (EQN.(19)) AND 

(C)MEDIATOR F m (EQN. (20)) COMPUTED FOR SOME FIXED 

VALUES OF THE PARAMETERS B =0.5, B =0.05, B =0.1 

^ ° r 

AND VARIOUS VALUES OF THIELE MODULUS ^ , tp] AND cp 2 m . 


(a) 



(b) 



(c) 



FIGURE 3: NORMALIZED CONCENTRATIONS OF (A) 
SUBSTRATES; (EQN. (18)) (B) OXYGEN F 0 (EQN.(19)) AND 
(C)MEDIATOR F m (EQN. (20)) COMPUTED FOR SOME FIXED 
VALUES OF THE PARAMETERS B 0 = 0.05 , B r = 0.1 AND $ 2 S = 1 
AND VARIOUS VALUES OF NORMALIZED SURFACE 
CONCENTRATION OF SUBSTRATE 5 S ,WHEN(I) B s = 0.001 (II) 

B s = 0.005 (III) S s =0.01 (IV) B s = 0.05 (V) S s =0.1AND 
(V) S =0.5 
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(a) 



(b) 



(c) 



FIGURE 4: NORMALIZED CONCENTRATIONS OF (A) 
SUBSTRATE^ (EQN. (18)) (B) OXYGEN F 0 (EQN.(19)) AND 

(C)MEDIATOR F m (EQN. (20)) COMPUTED FOR SOME FIXED 
VALUES OF THE PARAMETERS 5^ = 0.5 , AND tp] m 1 AND 

VARIOUS VALUES OF NORMALIZED SURFACE 
CONCENTRATION OF OXYGEN WHEN, (I) B 0 = 0.001 (II) 

B o = 0.05 (III) 5 o =0.1 (IV) B o = 0.5 (V) B o = 1 AND (VI) B o = 2. 


(a) 



(b) 



(c) 



FIGURE 5: NORMALIZED CONCENTRATIONS OF (A) 
SUBSTRATES; (EQN. (18)) (B) OXYGEN F 0 (EQN.(19)) AND 
(C)MEDIATOR F m (EQN. (20)) COMPUTED FOR SOME FIXED 

VALUES OF THE PARAMETERS B = 0.5 , B = 0.05 AND a 2 = 1 

s ° ^ m 

AND VARIOUS VALUES OF NORMALIZED SURFACE 
CONCENTRATION OF MEDIATOR(II) B m = 0.005 (III) B m = 0.01 
(IV) B m =0.05 (V) B m =0A (VI) B m = 0.2 AND (VI) B m = 0.5 . 
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FIGURE6: NORMALIZED CONCENTRATIONS OF SUBSTRATE 
F s (EQN. (18)) VERSUS THE NORMALIZED DISTANCE x , 
OXYGEN F a (EQN. (19))VERSUS THE NORMALIZED DISTANCE 
x WHEN B s = 0.5 , B 0 = 0.05 , B m = 0.1 AND cp] = 1 AND 
MEDIATOR F m (EQN.(20)) VERSUS THE NORMALIZED 
DISTANCE x WHEN B = 0.5 , B = 0.05 , B m = 0.1 AND (pi =1. 

S ° m 'rm 


(a) 



5 r = 0.2 



(c) 



(d) 



FIGURE7: VARIATION OF NORMALIZED CURRENT J ox WITH 

(A) NORMALIZED THIELE MODULUS FOR THE MEDIATOR (B) 
NORMALIZED SURFACE CONCENTRATION OF THE 
SUBSTRATE B s (C) NORMALIZED SURFACE CONCENTRATION 
OF OXYGEN B 0 AND (D) NORMALIZED SURFACE 
CONCENTRATION OF MEDIATOR B m . 

Conclusions 

The system of non-linear reaction diffusion equations 
of the model has been solved analytically. The 
approximate analytical expressions for the 
concentrations of substrate, oxygen, mediator and the 
current at the enzyme-membrane electrode geometry 
are obtained using the new Fiomotopy perturbation 
method. This analytical result is useful for 
improvingthe sensor design. The extension of the 
procedure to the model with reduced mediator in the 
bulk solution seems possible. These analytical 
expressions can be used for the optimization of the 
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thickness of the enzyme layer or Thiele modulus 
which produces significant change in both the 
magnitude of the current and the general behavior of 
the system. 
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Appendix A: Basic Concept of the Homotopy 
Perturbation Method 

To explain this method, let us consider the following 
function: 

D„(u) - f (r) = 0, reD (A.l) 

with the boundary conditions of 

f)u 

B 0 (u, — ) = 0, reT (A.2) 

on 

where D 0 is a general differential operator, B 0 is a 
boundary operator, / (r) is a known analytical 
function and T is the boundary of the domain Q . In 
general, the operator D 0 can be divided into a linear 
part L and a non-linear part N . The eqn. (A.l) can 
therefore be written as 

L(u) + N (u) - f (r) = 0 (A. 3) 

By the Homotopy technique, we construct a 
Homotopy v(r,p ) : Qx[0,l] — » 91 that satisfies 

H(v , p) = ( l- p)[L(v) - L(u 0 )] + p[D 0 (v) - f(r)\ = 0 (A.4) 

H(v, p ) = L(v) -L(u 0 ) + pL(u 0 ) + p[N(v) - f(r)] = 0 (A.5) 

where p e [0, 1] is an embedding parameter, and u 0 is 
an initial approximation of the eqn. (A.l) that satisfies 
the boundary conditions. From the eqns. (A.4) and 
(A.5), we have 

H(v, 0) = L(v) - L(u 0 ) = 0 (A.6) 

H (v, 1) = D 0 (v) - f(r ) = 0 (A. 7) 

When p= 0, the eqns. (A.4) and (A.5) become linear 
equations. When p =1, they become non-linear 
equations. The process of changing p from zero to 
unity is that of L(v)-L(u 0 ) = 0toD o (v)-f(r) = 0. We 
first use the embedding parameter p as a small 
parameter and assume that the solutions of the eqns. 
(A.4) and (A.5) can be written as a power series in p : 


This is the basic idea of the HPM. 


Appendix B: Analytical Solution of the 
System of Non-linear Equations (12)-(16) 
Using New Homotopy Perturbation 
Method [16-27] 

In this Appendix, we indicate how the eqns. (18)-(20) 
in this paper is derived. To find the solution of the 
eqns. (12) - (14), we construct the new Homotopy as 
follows: 


0-/0 


d 2 F s 
dx 2 


a 2 

(F 0 (\)B 0 + FJ\)B m ) F S {\)B S +F S {\)B S 
K +{F 0 (\)B 0+ F m (i)B m ) 

f M 


J(B.l) 


+ P 



v 


( 1-/0 



{F 0 Bo+F m B m )F a B a =Q 
(F 0 B 0+ F m B m )F s B s+ F s B s 
+(F 0 B 0+ F m B m ) )_ 

tf(B 0 +B m )B a 

(B 0 +B m )B s +B s +(B 0 +B m ) 


d 2 F s 
dx 2 

r \ 

(F 0 B 0 +F m B m )F s B s 
(F 0 B 0+ F m B m )F s B s+ F s B s 
y + { F oB 0 +F m B m ) , 



(B.2) 


dx 2 

f 


(1 ~P) 


do Km, Km, 


(K m„ + F m (i )B m ) f s (i)5, + k m s 

+(F 0 (\)B 0+ F m {\)B m ) 


(B.3) 


+ P 


d 2 F„ 

dx 1 


-~(Po 


FB FB 

0 0 s s 


(F 0 B 0 +F m B m )F s B s+ F s B s 

+(F 0 B 0+ F m B m ) 


= 0 
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pp pa a 

dx 2 (B 0 + B m ) B s + B s + (B 0 + B m ) 



PA PA 


(PA+ F m A)PA+PA 

+ { F o B o +F m B m) 


(1 -P) 


P PJVA PFA 

( F o m, + p m (i )B m ) f s (i )b s + f s m 

+ (F 0 (l)B 0 + F m (l)B m ) 



F B FB 

mm s s 


{PA+ F m B m )PA+PA 

+ ( F o B o + PAm) 




[ Pm <Pm A A 

dx 2 ( B 0 +B m ) B s +B s + ( B 0 + B m ) 



F m B m PA 


(F B +F B ) FB +FB 

V o o m m ) s s s s 

+ (F B +F B ) 

V V o o m m ) 

The analytical solution of the eqn.(B.2) is 

p =p„ +pp, +p 2 p, + 


Similarly the analytical solutions of the eqns. (B.4) and 
(B.6) are 

P=P +P F o+P 2 P 2+ (B.8) 


F m = F m„ +P F ml+P F m2 


Substituting the eqn.(B.7)-(B.9)into an eqn.(B.2) we get 


d 2 A + P F S +P 2 F s , +...) 


(P, +pP 2 +p 2 P, +AA 

A F m, + P F m, + P lF m 2 + -) B m 


\p +pP 2 +p 2 P, +AA 

v + (^m 0 + P F m, + P 2F m 2 + —) B m 

(P„ +pp, +p 2 P , +AA 


+(P„ +PP, +P 2 P , + -)A 


(P, +Pp 2 +P 2 p 3 + AA 
H F m n + P F m, +P 2 F m, + AA 


(B.10) 


Substituting the eqns. (B.7)-(B.9) into an eqn.(B.4) we get 


d 2 P + P F 0 , + p 2 p„ +...) 


( 1-/0 


PAP 

_ (A+A)A+A+(A+A)_ 

~d 2 (F + P F„ +p 2 F +...) 


(P 0 +pP, +p 2 P 2 +-)A 
(p +pp t pAp pPA 
\p +pP 2 +p 2 P, +AA 

y + ( F m„ + P F m, + P 2 P mi + -)B m 
(F +P F St + p 2 F +...)B S 


HP . +pF +p 2 F s , +...)B, 


(P, + P F 0 +p 2 F +...)B 0 


(B. 11) 


L l [Hp+pp+Pp+AA)) \ 

Substituting the eqns.(B.7)-(B.9) into an eqn.(B.6) we get 



~d 2 (F, t +P F Si +P 2 F Si + ...) 


'd 2 (F mo+P F mi+P 2 F m2+ ...) 

a-P) 

dx 2 

<p 2 (B 0+ BJB s 

(i -p) 

dx 2 

Vi B m B s 


_ (A + A)A + A+(A + A)_ 


_ (. B 0+ B m )B s +B s+ (B 0+ B m )_ 
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+Pl\ +P 2 J\ +-) 

dx 2 


+P 


-<Pm 


( F m c + P F m, + P F m 2 + -) B m 
\F 0 i +pF 02 +p 2 F 03 + ...)B 0 

+ ( F m„ +P F m t +P lF m 2 + -) B m , 

( F ,, +P F s, +P 2 f s 2 +-)B s 

+( F S „ + P F s, +P lF s 2 +-)B S 

F 0 , + P F o 2 +P 2 Fo, + •••)*„ 

+ (^m„ + P F m, + P lF m , + —) B m 


= 0 


(B.12) 


Comparing the coefficients of like powers of p in (B.10) 
we get 


P°: 


d 2 F 

^0 

dx 1 


-<Ps 


(B 0 +BJB s 


_(B 0 +B m )B s+ B s+ (B 0 +B m ) 
The initial approximations are as follows 

^ 0 (i) = i; F 'sS °) = o 

F s. (i) = 0; F\ (0) = 0 ,i = 1,2,3,... 


0 (B.13) 

(B.14) 

(B.15) 


Solving the eqn. (B.13) and using the boundary 
conditions eqns.(B.14)-(B.15), we obtain the following 
result: 


F 5 (x) = F So =^(x 2 -1) + 1 
Where Z:, 2 is defined in the text eqn. (21). 


(B.16) 


After putting the eqn. (B.16) into an eqn. (B.7), we 
obtain the solution in the text eqn.(18). 

Comparing the coefficients of like powers of p intoan 
eqn.(B.ll) we get 


P°: 


d 2 F 
^0 

dx 2 
-% 2 


B 0 B S 


_{B 0 +B m )B a +B a +(B 0 +B m ). 

The initial approximations are as follows 
F (1) = 1;F» (0) = 0 


(B.17) 


= 0 


(B.18) 


F 0 . (1) - 0; F' 0 (0) = 0 , / = 1,2,3,... (B.19) 


Solving the eqn. (B.17) and using the boundary 
conditions eqns. (B.18)-(B.19), we obtain the following 
result: 

F 0 (*) = F 0o J&- (x 2 - 1) + 1 (B.20) 

Where k \ is defined in the text eqn. (22). 


After putting the eqn.(B.20) into an eqn.(B.8), we 
obtain the solution in the text eqn.(19). 

Comparing the coefficients of like powers of p into an 
eqn.(B.12) we get 


d 2 F ' 


dx 2 


-<Pm 


« 

o B 0 +B m )B s+ B s +(B 0 +B m ) 


= 0 (B.21) 


The initial approximations are as follows: 

F m 0 (l) = 1 ;^ , m 0 (0) = 1 

B m . (1) = 0; F (0) = 0 


(B.22) 

(B.23) 


Solving the eqn. (B.21) and using the boundary 
conditions eqns.(B.22)-(B.23), we obtain the following 
result: 


F m (x) = F mo =^-(x 2 -x) + 1 (B.24) 

Where k 2 is defined in the text eqn. (23). 

After putting the eqn.(B.24) into an eqn.(B.9), we 
obtain the solution in the text eqn. (20). 

Appendix C: Scilab/Matlabprogram for the 
Numerical Solution of the Systems of 
Non-linear Differential Eqns. (12)-(16) 

function pdex4 
m = 0; 

x = linspace(0,l); 
t = linspace(0,100000); 

sol = pdepe(m,@pdex4pde / @pdex4ic,@pdex4bc,x,t); 
ul = sol(:,:,l); 
u2 = sol(:,:,2); 
u3 = sol(:,:,3); 

%plot(x,ul(end,:)) 

%xlabel('Distance x') 

%ylabel('ul(x,2)') 

%figure 
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%plot(x,u2(end,:)) 

%xlabel('Distance z') 

%ylabel('u2(x,2)') 

%figure 

%plot(x,u3(end,:)) 

%xlabel('Distance x') 

%ylabel('u3(x,2)') 

% 

function [c,f,s]=pdex4pde(x,t,u,DuDx) 
c = [1;1; 1]; 
f = [1; 1; 1] .* DuDx; 

Bs=0.5; B0=0.05; Br=0.1; 

al=sqrt(0.001); a2=sqrt(0.001); a3=sqrt(0.001); 

Fl=-al A 2/(l+(l/(u(2)^B0+u(3)^Br))+(l/(u(l)^Bs))); 


Appendix D: Nomenclature 


Symbol 

Meaning 

te] 

Total active enzyme concentration in the matrix 
(mmol/L) 

[EoA 

Enzyme concentration of the oxidized mediator 
(mmol/L) 

[£S] 

Enzyme concentration of the substrate (mmol/L) 

[ £ J 

Reduced enzyme concentration (mmol/L) 

[o,] 

Concentration of oxygen at any position in the enzyme 
layer (mmol/L) 

[0,1 

Oxygen concentration in the bulk electrolyte (mmol/L) 

[0,1 

Oxygen concentration in the bulk solution (mmol/L) 

P] 

Concentration of substrate at any position in the 
enzyme layer (mmol/L) 

PL 

Substrate concentration in the bulk electrolyte 
(mmol/L) 

[Med ox ] 

Concentration of oxidised mediator at any position in 
the enzyme layer(mmol/L) 

[Med ox i 

Oxidised mediator concentration in the bulk 
electrolyte (mmol/L) 


Oxidized mediator concentration in the bulk solution 
(mmol/L) 

D o 

Diffusion coefficient of oxygen (cm2 s-1 ) 

D s 

Diffusion coefficient of substrate (cm2 s-1 ) 

D m 

Diffusion coefficient of mediator (cm2 s-1) 

d 

Enzyme layer thickness (cm) 

y 

Distance from the electrode (cm) 

K,k 4 ,k 5 

Rate constants ( L mol-1 s-1 ) 

k 2 ,k 5 

Rate constants ( s-1 ) 

K o 

Partitioning coefficient for oxygen (none) 

K s 

Partitioning coefficient for substrate (none) 

K m 

Partitioning coefficient for mediator (none) 

B o 

Normalized surface concentration of oxygen 

B s 

Normalized surface concentration of the substrate 

B m 

Normalized surface concentration of mediator 

F 

Normalized surface concentration 

F o 

Normalized oxygen concentration 

F m 

Normalized mediator concentration 


F2 =- 

(a2 A 2^(u(2)^B0)/(u(2)^B0+u(3)^Br))^(l/(l+(l/(u(2)^B0+u( 
3)*Br))+(l/ (u(l)*Bs)))); 

F3 — 

(a3 A 2*(u(3)*Br)/(u(2)*B0+u(3)*Br))*(l/(l+(l/(u(2)*B0+u( 
3)*Br))+(l/ (u(l)*Bs)))); 
s = [FI; F2; F3]; 

% 

function uO = pdex4ic(x) 
uO = [1; 1;1]; 

O/ 

/o 

function [pl,ql,pr,qr] = pdex4bc(xl / ul,xr,ur / t) 
pi = [0; 0; ul(3)-l]; 

qi = [1; 1; 0]; 

pr = [ur(l)-l; ur(2)-l; ur(3)-l]; 
qr = [0; 0; 0]; 


Greek Symbols 


<P 2 s 

Thiele modulus for the substrate (Normalized) 

<Po 

Thiele modulus for the oxygen (Normalized) 

<p 2 m 

Thiele modulus for the mediator (Normalized) 


Subscripts 


o 

Oxygen 

s 

Substrate 

m 

Mediator 

ox 

Oxidized species 

red 

Reduced species 

t 

Total 

00 

Bulk solution 
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